Non-minimal coupling, exponential potentials and the w < —1 regime of dark energy 
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Recent observations and theoretical considerations have motivated the study of models for dark 
energy with equation of state characterized by a parameter w — p/p < —1. Such models, however, 
are usually believed to be inviable due to their instabilities against classical perturbations or poten- 
tially catastrophic vacuum decays. In this Brief Report, we show that a simple quintessential model 
with potential V{cj>) — Ae~'"^ and a gravitational coupling of the form (1 + |5|(^^)7l can exhibit, for 
large sets of initial conditions, asymptotic de Sitter behavior with w < —1 regimes. Nevertheless, 
the model is indeed stable at classical and quantum level. 
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I. INTRODUCTION 

The nature of the dark energy component responsible 
for the accelerated expansion of the universe is one 
of the most profound problem of Physics (for a recent re- 
view, see [!). The simplest way to describe dark matter 
is by means of a cosmological constant A, which acts on 
the Einstein equations as an isotropic and homogeneous 
source with equation of state — pa = /oa = A. Straightfor- 
ward questions about possible fluctuations of dark energy 
leads naturally to the introduction of a field </)(x) (the 
quintessence instead of the cosmological constant A. 
Besides the issue of fluctuations, the field description is 
preferable since for some modelsfthat ones, for instance, 
provided with tracker solutionsQ) the appearance of an 
accelerated expansion (de Sitter) phase is a generic dy- 
namical behavior, avoiding, consequently, problems with 
fine-tunning of initial conditions in the early universe. 
Quintessential equations of state are, in general, of the 
type prj, — wpiji, where the parameter w can vary with 
time. For minimally coupled models with usual kinetic 
terms, one has 



P0 (/)2 4- 2V{<1)) 



> -1. 



(1) 



According to the Einstein equations, cosmological models 
with accelerated expansion phases require source terms 
for which w < —1/3. 

Recent observational questions'^ and theoretical 
speculations^ have motivated the analysis of the pos- 
sibility of having realistic models for which, at least 
temporally, w < —1. Indeed, recent data from Hub- 
ble Space Telescope of Type la supernova at high red- 
shifts (z > 1)0 favorer scenarios with slowly evolving w 
and restrict its values to the range w = —1.02 (1q ^g)- 



They are compatible with a simple cosmological con- 
stant (w ~ —1), but if, however, one really has w < —1, 
both cosmological constant and minimally coupled scalar 
quintessence descriptions for dark energy are ruled out. 
Realistic models accommodating also the regime w < —1 
are, of course, welcome. 

In '91, Carroll, Hoffman and Trodden consider the vi- 
ability of constructing realistic dark energy models with 
w < —1. The so called phantom fields, i.e. minimally 
coupled scalar fields with "negative" kinetic energy, are 
usually invoked to construct such models. Despite that 
the dominant energy condition (in fact, all usual energy 
conditions) is violated in models with w < —1, the au- 
thors are able to introduce a classically stable model in- 
volving a phantom field. However, due to the peculiar 
kinetic energy of the phantom, the model is unstable 
against any quantum process involving it. In particu- 
lar, as the phantom potential is unbounded from below, 
there are catastrophic vacuum decays]^. Furthermore, 
some general results 10] suggest that any minimally cou- 
pled theory with w < —1 has spatial gradient instabilities 
that would be ruled out by CMB observations. These re- 
sults put severe doubts on the viability of dark energy 
models based on phantom fields. 

In this Brief Report, we notice that certain quintessen- 
tial models can exhibit a generic asymptotic de Sitter 
phase with many solutions for which w is slowly evolving 
and lesser than -1, without the introduction of any (clas- 
sical or quantum) instability. The quintessential field is 
assumed to be non-minimally coupled to gravity 
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-g{F{(f>)R-da4'd''(t>-2Vm, (2) 



where F(0) = 1 — (^cjp, £, < 0, and to have an exponential 
self-interaction potential 



(3) 
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Exponential potentials have been used recently in cos- 
mology, mainly in connection with tracker fields [Tll|. 
They appear naturally in higher-dimensional theories 
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and string inspired models Non-minimally coupled 
models are also commonly adopted. In 13], for instance, 
a quintessential model with ^ = 1/6 (conformal coupling) 
and V{(f>) = - was introduced. With such 

potential, however, conformally coupled models never 
exhibit asymptotic de Sitter phases. Moreover, mod- 
els with coupling F{(f>)R are generically singular on the 
hypersurfaces = 14], precluding the viability of 

quintessential models with conformal coupling. Some 
proposals to circumvent this singularity with the inclu- 
sion of higher ord er g ravitational terms in the action 
has been suggested!^, but the viability of the result- 
ing models is still unclear. In this work, we avoid this 
singularity by choosing hereafter ^ < 0. Non-minimal 
couphngs of the type = 1 -f — 0o) been 

considered in [T^ to study tracking behavior for poten- 
tials V{(t)) = M4+"/^, a > 0. The case of = 
was considered in n7\. As we will see, despite that 
the model given by l|2Jl is known to violate weak energy 
condition [13], the phase space for the comoving frame 
presents large regions of stability. 



II. THE MODEL 



where Fi{4>) = For the metric ©, the 

Klein Gordon equation reads 



G(0,0,g) 



+ Kff(0) = o, 



where 



and 



G(</), 0, H) = 3FiH + -(1 + 3F")F'^ 



{FV' - 2F'V) 



(10) 



(11) 



(12) 



The fixed points of the equations (|5|)- (|10|l are the con- 
stant solutions (j)(t) — (f) and H{t) = H, corresponding to 
de Sitter solutions a{t) cx e^* for which w = — 1. Despite 
that the potential |j2Jl has no equilibrium points, thanks 
to the non-minimal coupling, the model has indeed the 
fixed points {(j)-,±H^) and {(()+, ±H+), where 



-^(iTyTT^V4! 



(13) 



The Einstein equations obtained from the action Q 



are 



FG,b=Va</>Vfa</)-^(V,0V^ 



while the Klein-Gordon equation is 

D<j>-V' + ]-F'R = 0. 



2V ~ 2UF)+VaVbF, 
(4) 



(5) 



A relevant issue here is how to define an equation of state 
for the field 4> from the equations (@J and ©. The r.h.s of 
Q does not correspond to a covariantly conserved energy 
momentum tensor due to the presence of F[(j)) in its l.h.s. 
In order to define the pressure and energy for the field 
(/) in a consistent way with the continuity equation, one 
needs a covariantly conserved energy momentum tensor, 
and from Q we can get Gab = Tab with[l^ 



and 



Hl = 



3F (6+) 



(14) 



Tab = ^a4>^b 



-^(v,^v^ 



2y - 2nF)+VaVf,i^+(i- 

(6) 

Assuming an isotropic and spatially flat universe, 

d.'P = -dt^ + a^{t) {dx^ + dy^ + dz^) , (7) 

we get, from the temporal component of the Einstein 
equation the energy constraint 



The fixed points exist, of course, only for ^ < — tT^/4. The 
next Section is devoted to the study of the phase space 
of this model. As we will see, the fixed point (0_ , H-) is 
an attractor; large sets of solutions tend spontaneously 
to this de Sitter phase, irrespective of their initial condi- 
tions. When approaching the de Sitter point, solutions 
can have w < — 1, without any induced instability. 

We finish this Section with the definition of the pres- 
sure and energy . These quantities are defined from 
the energy momentum tensor as Tab = {p+p)uaUb+pgab, 
where Ua is a globally timelike vector. With the hypoth- 
esis of isotropy and homogeneity, we get, in the comoving 

frame, = 3iJ^ and p^ = — (2H + 3i?^^ , leading to 

-F)Gab- 



w 




(15) 



iH [FH + F'd)] = 



(8) 



As one can see, the ratio w defined in H15|) is not subject 
to the restriction w > —1. Note that the continuity equa- 
tion P(j, -\- 3H{pci, +Pci,) — 0, which is a direct consequence 
of the Bianchi identities for the Einstein equations, holds 
here. 



where H = d/a. From the spatial components, one gets 
the modified Friedmann equation 

-2FiH ^3{F + 2{F'f) H^ + —^ ,j)^-V- [H^ + V' 

(9) 



III. THE PHASE SPACE 



F' 



We study the phase space of the model by using the 
same semi-analytical approach used in [l3j ]. The phase 
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space is 3-dimensional {(j),(j),H), but due to the en- 
ergy constraint © the dynamics are restricted to a 2- 
dimensional submanifold. No restrictions are imposed on 
the {ip, (j)) phase portrait, but for the {(f), H) one, only the 
region > V/3Fi is dynamicaUy allowed (See Fig. 
Since H = is not allowed by the dynamics, the trajec- 
tories are confined to semi-spaces H > and H < 0. We 
are concerned here only with H > 0. We remind that real 
trajectories move on the 2-dimensional manifold defined 
by the energy constraint JH)). Therefore, in the projection 
on the plane ((/>, H) , each point on the allowed region 
corresponds, in fact, to two possible values for (p (two 
"sheets"), with the exception of the lines = V/3Fi, 
where only one value for (j) is allowed. Solutions aways 
"cross" from one sheet to another tangentially to the lines 
= V/3Fi. 
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FIG. 1: Tlie fixed points for the equations (151- 11011 on the 
plane {(j},H). The shadowed region (H^ < V/3Fi) is dynam- 
ically unaccessible. For this graphics, — ^ = a — 1. 

The aspect of K;ff(0) obtained from 112|l is crucial to 
the identification of the attractor points in the phase 
space. From (|12() . one has 



Kff (0) = Aa^e' 



1-^(1-6002 



(16) 



which can be integrated by parts in terms of the Ex- 
ponential Integral function i?i(x)[23|. From Fig. El it is 
clear that the fixed point corresponding to 4>+ is unstable, 
while (j)- should correspond to a stable one. However, 
despite the clear fact that 0_ is a minimum of Vcs , one 
cannot conclude safely that it indeed corresponds to a 
attractor due to the function G(0, 0, H) given by Hll|) . If 
G(0, 0, H) > 0, the solutions of the Klein-Gordon equa- 
tion around the point 0_ will be simple dumped oscilla- 
tions. This can be checked by introducing the function 





Veff 



Gicj),cj),H) = 3[F{cb)H 



FIG. 2: Aspect of the potential (the solid line) Vcs{(j)) for 
—1^ — a — 1. The contours correspond to the lines of constant 
L((j!), 0) 1171 1. Since they are closed around the fixed point (f>- , 
the condition L < along solutions implies that they tend 
asymptotically to (j>~. 



and noticing that L — —GcjP along the solutions of the 
Klein-Gordon equation (|l()|l . Provided that G{(f>, (p, H) > 
0, the function L{(j), (p) is a Lyapunov function [T^ for the 
fixed point 0_, assuring its stability. From the Eq. lO, 
we see that when 0(i) approaches H ^ and H{t) 
approaches establishing the attractor character of 
the fixed point (0_ , H_ ) . 
From Hll() . we have 

^—^F'{ct>)<pyi{F'm'H, 

(18) 

and from the energy constraint (jHJ, one has that F{(j))H+ 
F'{(j))(j) > on the semi-space H > 0, implying the posi- 
tivity of G and, hence, establishing the attractive charac- 
ter of the fixed point (0_,i/_), for ^ > —5/6. However, 
this is a very conservative lower bound for ^ Our exhaus- 
tive numerical simulations suggest that it can be consid- 
erably smaller. We could verify the attractive character 
of the fixed point {(j)-,H-) even for ^ < —100, suggest- 
ing that eventual amplifications forthcoming from the 
G{(j), (j),H) < regions are not enough to win the poten- 
tial Vcff around 0_ . A typical phase portrait is displayed 
in Fig. 13 corresponding to the case — ^ = A = a = 1. 
The attraction basin is considerably larger than the con- 
servative estimative base on the closed lines of constant 
L around the fixed point. Note that all solutions start- 
ing with Ip > <p^ are runaway solutions. Fig. 0| shows 
the curves w(i) for some solutions presented in Fig. |31 
According to (|15|1 . solutions approaching the fixed point 
from below have w < —1, while that ones doing from 
above have w > —1. 



IV. CONCLUSION 



+ Kff(0), 



We show here that non-minimally coupled quintessen- 

(17) 

tial models with exponential potentials can exhibit 
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FIG. 3: Typical phase portrait {(j>,H), corresponding to the 
case — 5 = A = a = 1. Due to the symmetries of Eq. lISl- llOL 
the H < trajectories are obtained from the H > ones by a 
reflexion on H — and reversion of the arrows (time reversal 
operation). 




asymptotic de Sitter phases for large sets of initial condi- 
tions. Some of these phases are characterized by a slowly 
evolving parameter w = Pcj,/P4, < —1, compatible, in 
principle, with the recent observational data. Our anal- 
ysis is based on the existence of a Lyapunov function for 
the Klein-Gordon equation that can be used to estimate 
the attraction basin of the relevant fixed points. As it was 
already mentioned, real attraction basins are typically 
much larger than these estimations. Our results are in 
agreement with the linearized analysis recently proposed 
in 13. 



We stress that the model presented here is free from 
the instabilities that are usually associated to phantom 
models. Classically, since F{(j)) ^ , the model is not 

tued with the anisotropic singularities described in 
Besides, since F{(^) is always positive and V{<^) 
is bounded from below, the model is also free from the 
quantum instabilities described in 9] for phantom fields. 

We conclude, therefore, that it is possible, in principle, 
to construct realistic models for dark energy with w < 
— 1. Relevant issues now are the introduction of matter 
fields and the study of the role played by the strong non- 
minimal coupling regime [13, I22| in the model. These 
points are under investigation. 



FIG. 4: The parameter w along some solutions presented in 
the phase portrait of Fig. According to l|15|l . solutions 
have w < —1 when H > 0, and ui > — 1 when H < 0. 
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